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Abstract. Our study is initiated by a multi-component particle system underlying 
the tiling of a half hexagon by three species of rhombi. In this particle system species 
j consists of lj/2\ particles which are interlaced with neigbouring species. The joint 
probability density function (PDF) for this particle system is obtained, and is shown 
in a suitable scaling limit to coincide with the joint eigenvalue PDF for the process 
formed by the successive minors of anti-symmetric GUE matrices, which in turn we 
compute from first principles. The correlations for this process are determinantal and 
we give an explicit formula for the corresponding correlation kernel in terms of Hermite 
polynomials. Scaling limits of the latter are computed, giving rise to the Airy kernel, 
extended Airy kernel and bead kernel at the soft edge and in the bulk, as well as a new 
kernel at the hard edge. 

1. Introduction 

A classical problem in random matrix theory is to compute the statistical properties of the 
eigenvalues, given the distribution of the entries of the matrix. For the well known Gaussian 
orthogonal, unitary and symplectic ensembles (GOE, GUE and GSE; see e.g. [For ) this problem 
can be solved exactly, with k x k determinant, and 2k x 2k Pfaffian formulas available for the 
fc-point correlations in the case of the GUE, and GOE, GSE respectively. Of these the GUE is 
special because determinantal expressions can also be derived for the correlations of the process 
specified by the eigenvalues of all the successive minors |JN06[[(7R06j . whereas for the GOE and 
GSE the correlations for the minors are not known. 

While interest in the eigenvalues of GUE matrices stems from its use in modelling the highly 
excited energy levels of classically chaotic quantum systems with no time reversal symmetry, the 
interest in the GUE minor process stems from a work of Baryshnikov [BarOlj relating to queues. 
Subsequently the GUE minor process has been identified by Johansson and Nordenstam within 
a number of random tiling problems [JN06| . In the language of stepped surfaces (which are 
equivalent to certain tilings) similar applications were noticed in ^ OR06j . while in |BFS07| the 
GUE minor process is encountered in the asymmetric exclusion process. 

One of the random tiling problems identified in |JN06| as relating to the GUE minor process 
is the tiling of a hexagon by three species of rhombi. This tiling problem is equivalent to non- 
intersecting lattice paths which at each step move up one unit or down one unit in a so called 
watermelon formation (see }KGVOO] ) . A very natural variation of the non-intersecting paths 
model is to impose the boundary condition that the paths cannot go below the x-axis |For89| 
IKGVO O. Gil03| , which in fact corresponds to the tiling of a half hexagon by the same three species 
of rhombi as referred to above (see Figure [1] below) . 

This variation initiates the study of the present work, where we begin in section 2 by analyzing 
such non-intersecting paths from the same viewpoint which in the study |JN06j gave rise to the 
GUE minor process. We obtain in corollary 12.111 a joint probability density function (PDF) 
which can be interpreted as a multi-component particle system in which the number of particles 
in species j consists of [j/2j particles (j — 1, . . . ,2n + I). The marginal distribution of species j 
can also be computed exactly (see theorem 12.81) , and is precisely that of a j x j anti-symmetric 
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GUE matrix (see e.g. [Fori Ch. 1]). By analogy with the findings from [JN06| . this motivates us 
to seek the joint eigenvalue PDF of the anti-symmetric GUE minor process, and to compare this 
against the joint PDF obtained from analyzing the non-intersecting lattice paths model. We do 
this in Section 3, and find the two joint PDFs are identical. 

It should be noted that the discrete processes that in this article and in [J NO 6] give rise to 
random matrix distributions is not the positions of the random walkers, but the holes between 
them, se section 2. It is also possible to obtain random matrix distributions by looking at the 
positions of random walkers directly, see for examble |Gil03llKMFW07] for examples of such 
models. For tilings of a hexagon, it is well known that the limit shape of the disordered region 
is a circle, |CLP98| . We believe that the limit shape in our half-hexagon model is a semicircle, 
but do not prove it. 

As already mentioned, the GUE minor process is determinantal, with the correlations between 
eigenvalues from prescribed minors being given by an r x r determinant. Explicitly, with the 
coordinate (s,y) denoting an eigenvalue from the s-th principal minor having value y, it was 
found that |JN06| 

P{r){.{{sj,yj)}j=i,...,r) = det[7^^u^'"((sj,yj),(sfc,yfc))]j-fc=i,...,^ 
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(here Hj{x) denotes the Hermite polynomials of degree j, satisfying J^Hi{x)Hj{ 
Sijnl2"y/Tr). In section 4, proposition 14. 3( we show that similarly for the anti-symmetric GUE 
minor process 
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In section 5 we turn our attention to various scaled limits of (jl.3p . In particular, two distinct 
soft edge scaling limits (neighbourhood of the largest eigenvalues of the respective minors), as 
well as a bulk and hard edge scaling (the latter referring to the eigenvalues in the neighbourhood 
of the origin) are computed. At the soft edge, with the species (i.e. minors) labels chosen to 
differ by 2n by a fixed amount, the limiting correlation kernel is found to be the well known Airy 
kernel (see e.g. jFori Ch. 4]) 



(1.5) 



K^°'\x,y) 



Ai{x + u) Ai{y + u) du 



independent of the species. The other soft edge scaling limit analyzed is when the species differ 
by 0{'n?/'^). In the case of the GUE minor process, study of this limit is motivated by a mapping 
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Figure 1. Tiling with rhombuses of "half a hexagon" and the corresponding 
non-intersecting random walks. The positions of blue particles are marked by 
triangles and the red particles by squares. The picture was generated with the 
method of [Pro03| . 



to the interface of a certain droplet polynuclear growth model |PS02j . Fluctuations of the latter 
are governed by the dynamical extension of the Airy kernel 

(16) K^°'H(r x)(r ^^ f C Ai(x + Ai(2/ + u) d., r, > 

(1.6) K j_jO^^-(..-..).Ai(,-,,)Ai(, + u)d., r,<r.. 

It is (|1.6p which again appears in this second soft edge scaling. In the bulk scaling of the GUE 
minor process, the limiting correlation kernel was found to be |FN08| 

(1.7) K^^^{{T,,x),{Ty,y))^{^{-) 

K ^ JR\[-1,1] 

This corresponds to the 7 = (isotropic) case of the correlation kernel for the bead pro- 
cess |Bou06j . It too is reclaimed in the appropriate bulk scaling limit of the anti-symmetric 
GUE minor process. The hard edge scaling limit has no analogue in the GUE minor process. 
The limiting correlations are given in proposition 15.51 

Note: While this work was being prepared, a research announcement by M. Defosseaux was 
posted on the arXiv, [DefOSj . stating results equivalent to our theorem 13.31 and proposition [331 

2. The half hexagon 

Consider p simple symmetric random walks started at (0, 2i — 2) conditioned to end up at 
(2A^, 2i — 2), for i = 1, . . . , p. They are conditioned never to intersect and never to go below the 
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X axis. Such configurations are in bijection to tilings of a half- hexagon, as seen in the Figure [TJ 
To count such configurations we need the following result of Krattenthaler et al. |KGVOO| . 

Proposition 2.1 (Theorem 6 in [KGVOO| ). Let ei < 62 < • ■ ■ < gp with Ci = m (mod 2), 
i ^ 1,2, . . . ,p. The number of stars with p branches, the i-th branch running from Ai — (0, 2i — 2) 
to Ei = (m, Ci), i = 1,2, . . . ,p, and never going below the x-axis equals 

A (e, + l)(m + 2z-2)! ^ 

Here, a star means a configuration of p simple symmetric random walks, or branches. Given 
that the red particles on line x = n occur at positions ei < e2 < • • • < Cp, the number of ways to 
tile the area to the left of line n is 

jj (e., + l)(n + 2z-2)! ^ ( (e^, + l^^ 

The number of ways to tile the area to the right of line n is 

^ (e, + l)(2jV-n + 2z-2)! yr ( f e, + 1^ f + 1^ 

Thus the probability of a certain configuration of red particles on line n is the product of the 
above two expressions divided by the total number of tilings of the area. 

Lemma 2.2. The probability distribution of red dots along line x = n is 

(2.4) V:Uei<e2<--- <ep] = 

Z-^ n ((e-+l)^-(e+l)^)^fT (e. + l)^(n + 2»^2)!(2iV-n + 2z-2)! 

ii<k ^^-+^^\n(i^+p)!(M^+p),(ii_^+p_i),(M^+p_i), 

for 0<ei<---<ep<2p + n~2 and Ci — n mod 2 for i — 1, . . . ,p. 

We are more interested in the distribution of the holes, i.e. the blue dots. This can be 
computed using a result from [Bor02j . To set this up, let X — {xi, . . . , xm} be a set of M real 
numbers. Consider measures on the following form. For w ; X — > R, 

m 

(2.5) = const J] 

l<k<l<m k=l 

is a probability measure which is supported on sets of to points. Notice that the measure above 
for the red particles is on this form. Denote the measure on the holes by 

(2.6) Vt'-"\A)^VZ\X\A), 

this is of course a measure supported on sets of M — m elements. Borodin proves the following 
formula. 

Proposition 2.3 (Proposition 2 in IBor02| ). Let u,v : X ^ M. be functions such that for all 

Xk e X, 

Then 

(2.8) = 
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With this tool we can now compute the distribution of blue particles. 
Theorem 2.4. On line x = n < N , for n even, the distribution of the blue particles is 

(2.9) nW{ei,e2,...,ep4 = 

np' ( 2N-n+e, , N|/ 2Ar-«-e, i „ _ i y 

,<,<,<,/^^ ^ l\{n + 2^-2y.{2N-n + 2^-m^+py.i^+p-ly. 

for p' = n/2 and < ei < 62 < • • • < Bp' < 2p + n — 2 and even for all i. For odd n it is 

(2.10) n\uc{ei,e2,...,ep4 = 

TT ffe - + 1)2-^6 +1)2)2 n 1^ 2 2 +P + 

ii<L' ^'^^'^^ fJ(n + 2^-2)!(2iV-n + 2^-2)!(^+p)!(H^+p-l)! 

/or p' = (n — 1) /2 and again 1 < ei < 62 < • ■ ■ < Bp' < 2p + n — 2 and 6^ orfd /or aZZ z. 
Proof. For n even, let X = {0, 2, . . . , 2Af — 2}. By proposition [231 we need to compute 



(2.11) n((:r + l)2-(y + l)2)-2 = 2- 



4M 



vex 



(M-f)!(M + l + f) '^ "2 
a; + 1 



For n odd, X ~ {1, 3, . . . , 2Af — 1} and we need to compute 

( {M - ^)\{M + ^^)\ 



(2.12) \{{{x + lf-{y + lf)-^^2-''' 



yex 



V (a; + 1)2 



Combining these with the formula in lemma [2721 gives the result. □ 

Let us fix some notation. For any k < n, let the positions of the [fc/2j blue particles on line 
k be denoted a;^'^'' > • • • > x''yki2\ ■ Let x^'^^ = {x'^i \ . . . , 2^[fe/2j )• For simplicity of notation later, 
let xp'^^-* = for ^ = 0, 1,. . . , L"-/2J- Obviously, for k odd, xf^ is always even and for k even, 
a;-'^^ is always odd. 

It so happens, as is easily seen from the picture, that the blue particles must fulfill certain 
interlacing requirements, 

(2.13) x^,<x^-^<x^ 

for all fc = 2, . . . , n and i = 1, . . . , [fc/2j . Henceforth we will write x'^'^'^'^^ >- x^^'^ when x^^^^^ and 
x^^'> satisfy the above interlacing. Let K{x"-) = {(x(^\ . . . ,a;(")) : x*^"^ >- a;*"^^^ y ■ ■ ■ y x'-^^. 

The next observation to make is that the distribution x^^\ . . . , x^"-*, given a;" is uniform in 
the cone given by the above inequalities. 

Proposition 2.5. The joint probability of all blue particles on the first n lines is 

^ ' 'biucia; ,...,x i- cardif(a;(")) ' ^ 

where Xfc(a;^'''\ a;('''+^)) = l{a;('^+^) >- a;^'^^} and cardii'(a;(")) is the cardinality of K{x'^^^). 
Seeing this, one would immediately want to compute cardiir(x'-"-'). 
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Lemma 2.6. Let Zi = 1, Z„ = Ilj^f n>i. For n even, 



(2.15) 

and for n odd (n > 3) 
(2.16) 



l<i<j<n/2 



cardK(a;(")) = Z,;! JJ 

l<i<j<(ri-l)/2 



Proof. The proof is an inductive one. With the products over i < j interpreted as unity in the 
cases n — 2 and 3, we see by inspection that both (12.15P and ()2.16p are correct in these base 
cases. 

Suppose (|2.15p has been estabhshed for n ~ 2N. Then we see that 



(2.17) 



cardi^(a;(2^+i)) = ^ cardi^(x(2^)) 

Use of the Vandermonde determinant identity in ()2.15p shows 



cardi^(x(2A')) ^ Z2"^det 



i,j = l,...,N 

The important feature is that each row in the determinant depends on a distinct coordinate, and 



(27V) f , , (2Ar+l) -, 

N+i-i from 1 to xXr+i_/ - 1 over 
odd values, or equivalently summing [x^^^l^i ~ l)/2 from zero to x^^'^^]^l/2 — 1 over integer 



allows the sum in (|2.17p to be carried out by summing each x 

odd values, ( 
values. Thus 



cari 



dK{x^^^^)^ Z-j^ dot 



(2N + 1) , 



^ {{x + i)xy-' 



x=0 



i,j = l,...,N 



The sum in column j is an odd polynomial in x'^-^^^^l of degree 2j — 1 with leading term 
(l/(2j — l)){xN+i-i/2)'^^^^ in row i. Adding appropriate multiples of columns 1, ... ,j — 1 to 
column j shows that this term can replace the sum, and so 

^ 1 

cardi^(x(2^^+i)) = Z-^ W ^ det [{x^^,^,/2f^-^]^^^^^^ 



3 = i 



.N ■ 



Removing a common factor of XN+i-il2 from each row and further use of the Vandermonde 
determinant identity gives (j2.16p . 

Analogous working shows that with (|2.16p established in the case n — 2N + 1, (|2.15p follows 
in the case n = 2N + 2. □ 

2.7. Asymptotics. 



Theorem 2.8. Under the resettling a;|"^ = -y2A^(l - 1/V3)a|"^ the measure V^^^^^ converges 
weakly to 



(2.18) 7'lSU{Al"\...,Ai";j = M/„"i n ((^. 



(")^2 _ ^\(»)^2^2rT.-(A<"')" 



l<i<j<n/2 



i/2 

He 

i=l 
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if n is even and 

(2.19) ^lGlE{Ai"\ . . . , Ag.j} = n ((4"')' - (^^' f)" 

l<j<j<Ln/2j »=1 

ifn is odd. Wn is the normalisation constant Wn — YilZi^^ Nn-21 where Nj is specified by |^.g7[ ) 
below. 

Proof. Apply Stirling's approximation to the formulas from theorem 12.41 □ 

As indicated by the use of the subscript aGUE, 'PaCUE happens to be the eigenvalue measure 
for the anti-symmetric GUE ensemble, as we shall see in the next section. We are interested in 
the full measure P^iuc^' and how its limit relates to aGUE matrices. 

Let AC^) = {X[''\. . . , a[^J,2j ) e (R+)", for fc = 1, 2, . . . . Also let Ap'+^^ = for Z = 1, 2, . . . . 

Consider the cone K(A(")) = {(A^^), . . . , A^")) : A^"' >- A("~i) > A^^)} where Af'^+i) A^*^) 

means that A^''^^^ > A^''^ > A^''^^^ > . . . . 

Theorem 2.9. Under the rescaling x'f '' = y^2iV(l - 1/\/3)a|"\ as N = p 00, the measure 
'^biue converges weakly to 

(oof)] ad) - xl(A''^A''^)■^^x«-l(A<"-^^A(")) („) 

l^-^Uj / aGUEmi^ , ■ ■ ■ , ' \ - VOlK(A(")) ' aGUEl-^ I 

where Xfc(A('=), A^^+i)) = IjAC^+i) >■ A^}. 

Proof. The convergence of a Riemann sum to an integral. The number of integer point in a cone 
K{x") suitably normalised converges to the volume of the cone K(A"). □ 

Lemma 2.10. Let Z„ be as in lemma W^ One has 

(2.21) volK(A(")) = Z-in((^i"^)' - i^ff) 

for n even and 

(2.22) volK(A(")) = Z-'W{{t^f - {\\-^f)J{{t^) 

i<j i 

for n odd. 

Proof. The proof of lemma 12.61 with summation replaced by integration. □ 
The interlacing conditions Xk can be written in terms of a determinant, 

(2.23) X.(AW, AC^'+i)) = dct[l{Af ^ < Af +^)}]i<,,,<M. 

(see e.g. |FR04[ lemma 1]) where Mk = \(k + 1)/2J and Ap'+^' = for Z = 0, 1, . . . . Also, 
introduce the standard notation for the Vandermonde determinant, 

(2.24) A(a^)= n 

l<i<j<n 

for a e M". This enables us to write the limiting measure on a very nice form. 
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Corollary 2.11. 

n-l njl 

(2.25) V^^^^^{\^^\ . . . , A(")} = 1^ A((A("))^) n det[l{Af < Af +^)}] J] e-^^'"'^^ 

" A:=l 1=1 

/or n even and 

(2.26) 7'lJi;;L{A«,...,A(")} = |^A((A(«))2)ndct[l{AW (A|"))e-(^^"')^ 

" fe=i 1=1 

/or n Olid. 

In the next section the measures T^aGUEm ^'"^ identified with the joint eigenvalue PDF for 
the minors of anti-symmetric GUE matrices. 



3. Eigenvalue pdf for the anti-symmetric GUE ensemble 

Consider the anti-symmetric GUE ensemble, which is the probability measure on purely imagi- 
nary Hermitian matrices with density Z~^e~"'"'^^ 1"^ . Here, Z is a normalisation constant. Equiv- 
alently, form a real Gaussian matrix with entries chosen independently from ^'^(O, I/a/S) and set 
H = ^{X — X^) . We seek to find the PDF of the eigenvalues of H and its principal minors. For 
this we adapt workings from |Bar01[|FR05j in which this task is carried out for H a GUE matrix. 
First, a technical lemma is required. 

Lemma 3.1. Let < ai < ■ ■ ■ < an be fixed real numbers. Let qi, qn be i.i.d. exp(l) 
random variables, specified by the PDF e^^ {x>0). Consider the random rational function 

(3.1) PW-^-t^- 

It has n positive zeros denoted < bi < ■ • ■ < bn, and their PDF is 

(") ft".'-'""' 

^ ' 1=1 

supported on ai < bi < 02 < ■ ■ ■ < On < bn- 

Proof. Claim 1: p has at least n positive roots. Since the qi are all non-negative, p(A) —^ —00 as 
A ai+ and p{X) ^ cxd as A — > a — . So there must be a root on each of the intervals (a^, 0^+1), 
for i = — 1. This also accounts for the inequalities oi < 61 < 02 < . . . < fen-i < cin- 

Further, since p(A) — > 00 as A ^ 00, there must be a root bn > On- 

Claim 2: p has at most n positive roots. The rational function p has simple poles at itoi for 
i = 1, . . . , n. Let 

n 

(3.3) Q(A)=n(A'-«D, 

1=1 

then pq is a polynomial of degree 2n + l that has the same zeros as p. The zeros already accounted 
for by the above argument are are ±&i, . . . , ±6„ and 0, so there can be no more. 
It follows from the above that it is possible to write 

(M p(A)-^nr..(A^-*f) 



nr.,(A^-<'?) 
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Comparing the residue at of p in (j3.ip and (13. 4p , an elementary computation gives that 



(3.5) 



nn / 2 



The PDF for the variables {gjf^i is 
(3.6) exp(-^(7i) 

and we want to change variables to {bi}f^i. The Jacobian J for that transformation is, up to a 
sign, 

(3.7) 



J = det 



6? 



Hi 



where the determinant is evaluated with the Cauchy double alternant identity. Inserting the 
expression for qi from (|3.5p simplifies this to 



(3.8) 

By expanding p.ip and 
(3.9) 



J = 2" 



Hi 



<i<j< 



at infinity and comparing the 1/A coefficient it follows that 

n n 
1=1 i=l 



Inserting this in p.6p and multiplying by the Jacobian (|3.8p gives the sought form 



□ 



Lemma 3.2. Let < ai < • • • < a„ 6e fixed real numbers. Let qi, g„ be i.i.d. exp(l) 
distributed random variables and let qo be r(l/2,l) distributed, having PDF (7rx)~^^^e~^ for 
a; > 0. Consider the random rational function 



(3.10) 



A 



E 



Mi 



A2 



It has n + I positive zeros denoted < bo < ■ ■ ■ < bn and their PDF is 



(3.11) 



2"+i A(b^) 
A(a2) 



j=0 



supported on < bo < ai < bi < ■ ■ ■ < On < bn- 
Proof. It is convenient to introduce oq = 0. Instead of 



choose 



(3.12) 



q{X) = Xl[{X'~a^). 



The the proof of lemma 13.11 goes through virtually unchanged with indices starting from zero 
instead of from one. The Jacobian expression from (j3.8p can then be simplified as 

/oio\ T _ on+l llo<i<j<nibt ~ ^j) -pT ^ _ f i ^« llo<i<j<nibi ~ ^j) Ili^O ^» 

1^0. i^j J~z 1-^2 _ 2^)11"'-^ t- TT ^o2_„2^^-r« 2 

ll0<i<i<nV"i "'j' i=o I I0<i<j<ny"'i "jVlli=l"i 

Computing the residue at the origin of p in (j3.10p and p.4p gives 

n n 

(3.14) n^'^-^oii"'- 



i=0 
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The expression corresponding to p.9p is 

n n n 

(3.15) ~E* = E«'-E^^ 

i=0 i=l 4=0 

The PDF for the variables {9i}"^o 

(3.16) !^P(-Er=o'?0 



V'^'Jo 

Muhiplying this with the Jacobian (|3.8p , inserting (|3.14p and (|3.15p gives the sought form p.lOp . 

□ 

Theorem 3.3. Let H be an n x n matrix from the anti- symmetric GUE ensemble. Let Hk 
be the k x k principal minor of H. Let A'-'"'^ = (aJ'^\ . . . , Aj^^^^gj) positive eigenvalues of 

Hk, ordered so that A*^ > X^^i- Then the joint PDF of A^^-*, . . . , A*-"-* is precisely that given in 
corollary \2.11\ 

Proof. Such a matrix H has the property that if A is an eigenvalue of H , then so is —A. Also, if 
the size of H is odd, this implies that one eigenvalue will be zero. 

The proof is an inductive one. A 2 x 2 matrix from this ensemble is of the form ( _?q g ) where 
a G N{fd, l/\/2). Its eigenvalues are ±a, confirming the theorem in the case n — 2. 

First, let n be even. Consider an n x n matrix A from this ensemble. The induction assumption 
is that its eigenvalue PDF is known. Consider the (n+ 1) x (n+ 1) matrix given by bordering A, 



(3.17) 



A w 
w* 



Here, w is a column vector of n purely imaginary numbers, all N (0,1/^/2). The star means 
transpose and complex conjugate. 

The eigenvectors of A can be paired up in the following way. If v is an eigenvector cor- 
responding to eigenvalue A then v is an eigenvector corresponding to eigenvalue —A. Con- 
sider a normalised eigenvector, \v\ = 1. Since v and v must be orthogonal to each other, 
I Rewp = j{v + v,v + v) — ^, where (•, •) denotes the inner product. 

Let C = [vi,vi,V2, ■ ■ ■] be the matrix whose columns are all the eigenvectors of A. Then 



(3.18) 



C* 0\ f A w\ fC 0\ _ f D C*w 
1 [w* [o Ij [w*C 



where D is a diagonal matrix of the eigenvalues. It follows from the above considerations of 
eigenvectors and an elementary calculation that ■w*C = (ai,ai,a2,...) where each a; is a complex 
number, the real and imaginary part of which are N{0, l/\/2). Let p„(A) be the characteristic 
polynomial of A and say that the eigenvalues of A are ±/xi, . . . , ±/i„/2- Of course the eigenvalues 
of A give the factorisation of pn as 

(3.19) Pr.{X)^{fil-X^)...{iil/,-X'). 

Then, it can be shown, say by expanding along the last row of the RHS of p.lSp . that the 
characteristic polynomial of that larger matrix is such that 

n/2 

(3-20) ^^ = A-y:^^ 



A computation shows that 20^0^ is exp(l) distributed. So we now need to find the PDF of the 
zeros of this random rational function, which is precisely what is given by lemma [XTl Multiplying 
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the expression that the induction assumption gives us for n with the conditional PDF from 
lemma [XT] proves the statement for -ffn+i when n is even. 

Assume now that n is odd. Do the same construction but the matrix A will now have one 
eigenvalue which is zero. Performing the same bordering as in (|3.18p . only this time w*C = 
(oi, ai, . . . , a„, a„, jfe). where h is Af(0,l/2). As above the characteristic polynomials for the 
n X n-matrix and the (?i + 1) x (n + l)-matrix are related by 



(3.21) 

Apply this time lemma 



(«-l)/2 



2a, a," 



Pn+i(A) ^ 

to prove the statement for Hn+i for n odd. 
4. Correlation functions 



□ 



4.1. The method of Borodin and collaborators. Consider 2n + 1 discretisations A^i, . . . , 
A^2n+i of the half interval [0, oo) each containing the point 0. On each set Mj distribute points 

,^L2j + l/2j „,;^u „2;-i 



W^\t!i^'^^ with x: 

signed) measure by 



= 0. On the configuration of points X := Uj=T^{^i} define a (possibly 



1 " 

(4.1) -Wdei [W2i-i{ 



xf-\xf) 



det \W2i{xl',x- 



1=1 



xdet fe-i(4"+i)] 



for some functions {Wi : Mi x Mi+i '^}i=i,...,2m and {qi : M2n+i ^}i=i,...,n- After 

making use of the Vandermonde determinant expansion one sees that the measures in Corol- 
lary [2TTl] are of this general form. One viewpoint is that (14. 1|) specifies a multicomponent system 
in which species Z (? = 1, . . . , 2n + 1) consists of [(/ + 1)/2J particles. But for / odd one of these 
particles is fixed at the origin, so the number of mobile particles for species I is [//2J- We seek 
an explicit formula for the correlations of the general multi-component system specified by ()4.ip . 

In fact a generalisation of this very problem has been addressed and solved in a recent work 
of Borodin and Ferrari, |BF07| . This work extends earlier work of Borodin and collaborators 
|BR05[[EFPS07j . However all these developments are relatively new, so it will be to the benefit 
of the reader that we give an account of the derivation herein rather than just quote the result 
(this is further justified as the proof of theorem 4.2 in [BF07| is very brief and calls for a detailed 
knowledge of the relevant results from [BROSllEFPSOT] 

To begin introduce the |A^2n+i| x n matrix = [qj-i{xi)]xii^M2n+i,3=i.---,n], the |A^2/-i| x 
|A^2; I matrices W2;-i [W2i-i{xj,yk)]x,eM2i-i,ykeM^i and the |7W2z|x |7W2;-hi I matrices 1^2/ := 
\W2i(xj,yk)\xj^M2i,yk<^M2i+i- Further introduce the n x \Ai2i\ matrices Ei with entries in row I 



all ones and entries in all other rows all zeros. Finally with M 
define the \M\ x \M\ matrix L by 



{l,...,n}UMiU---UM2n+i 



(4.2) 









El 





E2 .. 


En 











-Wi 























-W2 














































. -W2n-1 























-W2 























for block zero matrices of appropriate dimension. 
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Note that the rows and columns of L are labelled by the elements of the set AA in order. For 
Y a subset of M , introduce the notation Ly to denote the restriction of L to the corresponding 
rows and columns. Then, in accordance with the general theory of L-ensembles (see |BR05j ) one 
sees that the measure (14.11) can be written 



, . detL{i,...,„}uAr 
^ ' det(l* + L) 

where 1* is the A4 x Ai identity matrix with the first n ones set to zero. The significance of this 
structure is the general fact that the correlation function for particles at F G A4i U • • • U A^2n+i 
is given by [BR05j . 

(4.4) p{Y) = detKY, K=l-{l + LM\{i,....r})-^ 

for 1 the identity of appropriate dimension. The correlation functions are thus given by a 
determinant of size \Y\ — the number of particles specified in the correlation — and so by definition 
the measure (|4.ip is a determinantal point process. 

Following [BFPS07llBF07| we will now proceed to isolate special structures in the functional 
form (j4.ip which when present allow (1 + L)~^ in ()4.3|) to be computed explicitly. For this 
purpose, introduce the matrix 



(4-5) W^K,) = 



W,...Wj.i, i<j 

For p = 1 , . . . , 2n and j = 0, . . . , n — 1 set 

(4.6) K-2j-2ix) = (M^[p,2„+i)*).,„-j, x&Mj> 

(we use the notation (^)a,fc to denote the element in row a column b of the matrix A). Define 
functions {*i'p_2/-2} '=o,i constructed from 

(4.7) Span{{(£;iT4^[2»,p))2i,a:}»=l,...,[(p-l)/2] U {(5p,ovon£'p/2}} 

by the orthogonality requirement 

xeMp 

for j, fc = 0, 1, . . . with p - 2j - 2, p - 2k - 2 > 0. Also set ■^f'+'^{x) = qj{x) (j = 0, . . . , n - 1). 
In terms of this notation, theorem 4.2 in |BF07| as applied to the measure ()4.1|) isolates on a 
special structure which when satisfied allows the elements of Ky in (|4.3p to be made explicit. 

Proposition 4.2 (Theorem 4.2 in |BF07j ). Define the [m/2] x [m/2] matrix Bm by 

[m/2] 

(4.9) {E^W[2^,„^))2^,x + <5^,m/2 (£^m/2 )m,x = i^^n)t,l'^>'^~2l{x)- 

1=1 

Suppose Bm is upper triangular. With the set Y in J^.gp labelled {(s;,y;)} where yi G Msi we 
then have 

ls/2] 

(4.10) Kiir,x),{s,y)) (ify)(.,.),(s,y) = -(VK[,,,)),,y + ^ *;_2i(a;)*s-2;(y)- 

1=1 

Proof. We follow both the proofs of jBFPS07[ lemma 3.4] and [BF07[ theorem 4.2], and make 
use of jBROSj . The task is to compute the inverse in 
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Write L in the structured form 

(4-11) ^=[c D-1 

where 

(4.12) B=[0,Ei,0,E2,...,0,En,0] 

(4.13) C=[0,...,0,vl.]^. 
We know from [BROS) lemma 1.5] that 

(4.14) K = 1- D~'^ + D~^CM-^BD-^ 
with M := BD^^C and furthermore 

(4.15) D-^ = l + [VF,,,)]„=i,...,2„+i. 
From this latter formula we compute 

(4.16) D-'C = [iy[l,2„+l)«', . . . , l^[2n.2n+l)«', ^-J^, 

and we compute that the m-th member of the block row vector BD^^ is equal to 

[m/2]-l 

(4-17) ^ EkW[2k,m) + Srn,cvcnEjn/2- 

k=l 

These formulas in turn tell us that the (j, m) block of K is equal to 

[m/2]-l 

(4.18) - +W^[j,2n+1)*M~1( ^^tW^pfe^m) +'5m,cvcn£^m/2), 

k=l 

where W[j^2n+i) is to be replaced by 1 for j = 2n + 1. They tell us too that 

ri-l 

(4.19) M= (^i?,VF[2fe,™))*. 
Define the (2n + 1) x [to/2] matrix ^™ specified by 

(4.20) (ci>™).,.^"-(-")' l^^^[-/2] 

I 0, otherwise. 

The orthogonality (|4.9p used in conjunction with (|4.18p shows 

[m/2]-l 

(4.21) ^ -EfeVFpfe^m) + (5„i,cvcn-£'m/2 = 

fc=l 

and furthermore M = i?2n+i- These formulas, together with the assumption that B„i is upper 
triangular, give that (|4.18p reduces to 

(4.22) - W^[,-™) +W^[,-2„+i)«'*™- 

It follows from this and the definition (|4.6p that the element in row (r, a;) column (s, y) is equal 
to (|4.10p . as required. □ 



Bm 




As formulated proposition 14.21 applies in the setting that the domains have been discretised. 
Minor modification of the definitions gives that (|4.10p remains valid in the continuum limit. For 
this purpose, introduce the notation 



(4.23) {a b){x,y) — / a{x, z)b{z,y) dz 



(4.24) W^,,,){x,v) 



14 PETER J. FORRESTER AND ERIC NORDENSTAM 

and note that such an operation is the continuum Umit of matrix multiphcation. Then specify 

\Wi* ■ ■ ■ *Wj-i){x,y), i<j 
0, i>j 

(4.25) *p-2,-2(a^) = iW[p^2n+l) * qn-j-l)ix) 

and for p = I, . . . ,2n define fmictions {^p_2i-2} i=o.i,.... constructed from 

(4.26) SpanJj/ W[2^,p))it, x)\ U {5p,evcn} i 

I J »=l,...,[(p-l)/2] J 

by the orthogonahty requirement 

poo 

(4.27) / 'f;^2j-2ix)^l-2k-2i^)dx^S,,k. 

Jo 

Consider now the exphcit form of these quantities for the PDF in coroUary 12.111 There, 
independent of i, Wi{x,y) — Xx<y and so 

(4.28) W[,,,^{x,y) = ^—j-^x.<y{y " a^)'"'~'. 
Furthermore, from the Vandermonde determinant identity 

n 

(4.29) Y[x, Y[ (4-^) =det[a;|'=-i]j-fc=i,...,„cxdet[p2fc-i(a;j)]j,fc=i,...,« 

4=1 l<j<fc<n 

vahd for any polynomials pi{x) of degree I. We choose the pi{x) so that 

(4.30) / e'"" P2j-i{x)p2k~i{x)dx (X Sj^k- 



This is achieved by setting pi{x) — Hi(x), where Hi(x) is the Hermite polynomial. Thus qi{x) 

2 

e^^ H2i^i{x) and so 

(4.31) ^^2,-2(^)= e-^'il2„-2,-i(y)(y-x)2"-J'd2/. 
Making use of the Rodrigues formula 

(4.32) H,{x) = {-iye^'-^e-^' 
it follows from this that for j > 

(4.33) ^P^ix)^e-^'H,ix) 
while for j < 

(4.34) = {y^x)-^-'e-y'dy. 

It follows from (|4.28p that (|4.26p has the explicit form 

(4.35) span{{a;P"2'},=i^...^[(p_i)/2] U {5p,cvcn}} 
From this set we can construct 

(4.36) '^l-2^i^) = J^Hp^2;{x) 
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for J = 1, 2, ... , where 

/•OC 

(4.37) {Hj{x)fe-''^dx^V^2^-^j\ 

Jo 

which indeed exhibits the orthogonahty (j4.27p . Furthermore, the continuum hmit of the change 
of basis formula (|4.9p reads 

[m/2] 

(4.38) + <5,™/2 = E (5m).,;*™-2/(2^) 

1=1 

and this with ^m-2ii^) ^ polynomial of degree m — 2Z as in (|4.36p gives that B,n is an upper 
triangular matrix. We are therefore justified in appling proposition 14.21 

Proposition 4.3. Let '^^{x) be specified by and |7]^p for j > and j < respectively. 

Let ^^{x) = -^Hj{x) so as to be consistent with Ii4.36\ ). The r point correlation is given by 

(4-39) P(r)((sj,2/j)j=i,...,r) = det[K{{sj,yj), {sk,yk))]3,k=i,...,r 

with 

^ [*/2] 

(4.40) i^((s,x),(t,y)) = -— — — x.<,(y-x)*-^-^ + E*^_2;(^)$U(y) 

^ ' 1=1 

Equivalently, for s > t 

(4.41) i.((.,.),(t,,)).e-^^E^^^)^^^ 

1=1 ^ 

while for s < t 

(4.42) Kiis,x)At.y))^-e-^' E ^^^'^^^^^-^'^^^ 

/ — — CXD 

Proof. The only remaining task is to derive (|4.42p . For s < t consider 

(4.43) - ^^_l_^y X.<\y\{sgnyY{\y\ - xf-^-^ 

which for ?/ > is equal to the first term in ()4.40|) . For t even (odd) this is an even (odd) 
function of y and so can be expanded in the basis {H2j+e{y)}j=a.i,... where e = 0, 1 for t even, 
odd, according to 

(4.44) ^ ' f;^g^±iM re-\u-xr^-'H,,^Au)du, 

valid for x > 0. Making use of the Rodrigues formula (|4.32|) and recalling (|4.34|) this can be 
rewritten 

(4.45) -e-^' E Hs-MH,^My) _ j2 ^U,i^.)<fUM 

fc=-oo * p=[s/2] + l 

Substituting this for the first term in (|4.40p and writing in the resulting expression 

ls/2] ls/2] 

(4.46) E ^U,,ix)^Uiy) = e~^' ^ ^-^^■(;)^--2^(^) 
1=1 k=i 

gives (0321). □ 
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4.4. The method of Nagao and Forrester. In the recent work [FN08) . the correlations for a 
class of multicomponent systems generalizing the GUE minor process, in which species / consisted 
of I particles constrained so that they interlace with the particles of species l + l, were computed in 
two ways. One was by using the generalized formula from BFPS07, lemma 3.4], while the other 
adapted an approach to multi-component determinantal processes due to Nagao and Forrester 
|NF98| . In this section we will show how proposition 14.31 can be reclaimed by making use of this 
latter method. 

With W,{x,y)^W{x,y) 

— Xx<y, the first step is to rewrite (|4.ip so that it reads 



(4.47) 



C 



Y[det 



1=1 



X det 



L(n-i)x(n-i) 











'-(n-l)x(n-l) 



j=l,...,l 

[l]fe=l,...,i 







{n — l)xl 



0,x(„-o [W{xf,xf+%,j=i,...,i 
xdetfe-_i(4"+i)],-fe=i,..., 



Here ki{x) is arbitrary as the determinant does not depend on k;, and C is a normalization which 
may vary from equation to equation below. Furthermore, qj^i{x) = H2j-i{x). Proceeding 
as in the derivation of the equality between (|4.43p and (|4.45p we deduce 



(4.48) w{x,y)^— e-'^dt + e-^'Y. 



H2k-l{x)H2k{y) 



N. 



2k 



or alternatively 
(4.49) 



W{x,y) 



E 

fc=i 



H2k{x)H2k+l{y) 



2k+l 



both valid for x,y > 0. We substitute the first of these in (|4.47p . after choosing 



(4.50) 



dt 



therein, and the second of these in the second of the determinants in (|4.47p . 
Introduce now the notation 



(4.51) 

so that {ri2j{x)}j=o,i, 
Set 



-H,ix), jr-^N. 



1/2 



and {?72j+i(a;)}j=o,i,... each form a set of orthonormal functions on [0, oo). 



(4.52) 
(4.53) 



oo 



72fc 
72fe 



fc=0 



72fc+l 



V2k-i{x)r]2k{y) 
V2k{x)r]2k+i{y)- 
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From the above working we then see that (I4.47P can be written 



(4.54) 



L(n-Ox(n-0 







X det 



lx{n—l) 



L(n-i)x(«-0 



[l]fe=l,....i 

[0(°)(xf-\xf)],.i,...,_i 

j=l,...,l J 







(n—l) X I 



xdet[iij^i{xl"+%^k=i,...,n 



To proceed further, set 
(4.55) 

and use this to define 
(4.56) 



(,)_ J77,(a;[^)), J>0,;>1 
1 Sj,2i-i, otherwise 



A 



(4.57) 



G 



(s,t) 

J,; 

s.t) 
3,1 



E l2k+s (s) (t) 
, ^ '^2fe+sj-n+[s/2]'?2fe+t,;-n+[t/2] 

oo 

E 



72fc+i 

72fc+a (s) (t) 
^2;,+j'^2fc+sj-ri+[s/2]'?2fc+t,;-ri+[t/2]- 



One then sees that 

(4.58) det[r;2,_i(a:f"+i)],-z=i,...,i (x det[Af;+i'i)],- z^i,...,„ =: det A^^^+LD 



(4.59) 



det 



[l]fe=l,...,Z 



J=l,...,i 



cx det[G 



(2/-l,20i 



\j,k=l,...,n 



detG(2'-i^2/) 



(4.60) det[</.(^)(xf ,a;f+i)],,=i,..., cx det[Gj^^'^'+^)],,fc.i,...,„ =: det G(2'^2'+i) . 

Consequently it is possible to rewrite (|4.54p as 



(4.61) 



Idet 



■^(2n+l,l) ^(2n+l,2) ^(2ri+l,3) 

-G(i'2) 

-G(2.3) 















^(2n+l,2ri+l)' 
_Q{l,2n+l} 
_^(2,2n+l) 

_^(2n,2n+l) 



This has the same structure as [FNOS) eq. 4.23]. Simple manipulation detailed in this latter 
reference shows that this in turn can be rewritten to read 



1 



-det[F(^'*)],,«=i,...,2„+i, = 



A^''*\ s>t 
B'-'^-^K s<t 



(4.62) 

where B'-"-*^ := A^"'*) - G'-'-*\ Noting from the definitions of yl'"^*) and G^'-*^ that for s > t, 

reduces to 
(4.63) 



f]"'*^ cx 5j^i for j,l <n - [s/2], while for s < t, F^f'' ^ 0, j < n - [s/2] ov I < n - [t/2], this 



J;det[/(^'*)],,t=i,...,2„+i 



G 
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where is the s x t-matrix with entries 

^ " ^i-'[i/2]+n,i-[t/2]+n 

(4.65) 

These exhibit the reproducing property 



It-Ik 

't-2k 



(4.66) ftr'f^'d^ 



^(s,t) fit.u) ^Jt) _ / /j,™"^ ' S>t>UOT S<t<U 



0, otherwise 



which aUow the integrations required to reduce (|4.63p down to the r-point correlation function 
to be performed (see e.g. |For j ) . giving 

(4.67) P(r)i{{Sj,X,)}j = l.,...,r) = det[f^^^''''\,k=i,...,r 
Substituting (|43T|) in (|465l) shows 

(4.68) f^f ^e-^'/'K{{s,x,),(t,xi)) 
and this the result (|4.39p of proposition 14.31 is reclaimed. 



5. Scaling limits of the correlations 

In the recent study [FN08| various scaling limits of the GUE minor process were computed. 
Two involved species which differed by a fixed amount. In this situation the eigenvalue coordi- 
nates were chosen to correspond to the neighbourhood of the spectrum edge in one case, and the 
bulk of the spectrum in the other. A further scaling limit, in which the species differ by 0{N'^^^) 
and the eigenvalue coordinates correspond to the spectrum edge, was also analyzed. 

In the first of these situations, the scaled correlation kernel was computed to be the Airy 
kernel (jl.Sp independent of the particle species. In the second situation the correlation kernel 
for the isotropic bead process ()1.7|) was found. In the final situation, the dynamical extension of 
the Airy kernel (|1.6p was obtained as the scaled correlation kernel. 

The anti-symmetric GUE minor process of the present work permits the above three scalings. 
We will show below that at the soft edge the limiting correlation kernels (|3.5p and (|1.6p are 
reclaimed in the cases of the species differing by 0(1) and by 0{n'^^'^) respectively. In the bulk 
it will be show that the limiting correlation kernel l|1.7p is reclaimed. The antisymmetric GUE 
minor process also permits the well known hard edge scaling |For93j , involving eigenvalues in the 
neighbourhood of the origin. For eigenvalues from the same species, which is an even label, the 
limiting correlation kernel is 

(5.1) K'^{xTy):—2 / sin(7rwx) sin(7ru?/) 

Jo 

while for odd labelled species it is 

(5.2) K^{x,y):=2 / cos^nux) cos{nuy) du. 

Jo 

In the case of hard edge scaling of the correlation kernel with variable species label, differing by 
a finite amount, a generalization of the integral representations (|5.ip and (|5.2p is found. 
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5.1. Soft edge scaling. In the N x N GUE the soft edge scahng corresponds to the change of 
eigenvalue coordinates 

(5.3) y^ = V2N- ^' 



which has the effect of moving the origin to the neighbourhood of the largest eigenvalue and 
scaling the distances so the inter-eigenvalue spacings in ths neighbourhood are order unity. The 
same soft edge scaling (|5.3p applies for (2n + 1) x (2n + 1) anti-symmetric GUE matrices except 
that N h-^ 2n. We consider first this scaling with the species differing from 2n + 1 by a constant, 

(5.4) Si = 271 + 1 - a 

Proposition 5.2. For the soft edge scaling specified by I15.3\) with N i—^ 2n and |5.^[), 



(5.5) _^^-,„„„,),(..»))^^A-m,«) 

where K'^°^^ is given by iS. 5]) and a„(c, = e"''^"-' ' ^(2n)^'^/^. Hence we have the pointwise 
limit 

(5.6) lim f-=i— ) p{{{sj,y,)}j=i,...,r)=dei[K^°'\Y,,Yk)],,u=i,...,r 



n—yoo 



independent of the particle species. 

Proof. Substituting (|4.36p in (|4.40p shows that for Sj > s/ (cj < q) 

2e-y^ 1 

(5.7) K{{sj,yj),{si,yi)) ^ ^^i-afcfg _ 2^)1 ■^^^-^2fc(yj)gs,-2fc(2/i)- 

k—i 

As done for the derivation of the corresponding result in the case of the GUE, |FN08| . we make 
use of the uniform large N expansion |01v74j 

(5.8) e'^'/^HM-kix) = 7ri/42(^-'=)/2+i/4((iv - k)\y/^N-'/'^x 

k \ , ^..,-2/3^ f0(e-'=/^'''), k>0\ 



0(1), fc<0 



in the summand of (|5.7p and so obtain 

(5.9) K{{s„y,),{si,m)) ~ 2e-(2»)^^^(^.-^'^)2-(^^-')/2V2(2n)-V6x 



^/(2n+l-c,-2A:)!V/' 2k 2k . 



fe=i 



We observe that the leading contribution to the sum comes from terms k = 0{n^^'^). In this 
regime 



(2n + 1 - Cj - 2ky. 



1/2 



(5-10) 77;— — Trrir ~ (2?^) 

^ ' \{2n+l-ci-2ky.J ^ ' 

and so this term can be factored out of the summand, leaving us with a Riemann sum approxi- 
mation to the definite integral (|3.5p which implies ()5.5p 

We know from |FN08| that in the case Sj < si [cj > q) the form (14.42p is not appropriate, due 
to the resulting Riemann sum not being convergent. Instead, we make use of (|4.40p . breaking the 
sum over I therein into the ranges / £ [1, [s/2]], I e [[s/2] + 1, [t/2]] to deduce that for the given 
scaling K((sj,yj), (s/,?//)) is to leading order given by the right hand side of (|5.7D . but with the 
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upper terminal replaced by [s/c/2]. This latter difference does not effect the leading asymptotic 
form, and so (|5.5p is valid in all cases. □ 

We turn our attention now to soft edge scaling with the species separated by 0((2ri)^/^). 

Proposition 5.3. Scale Si according to 

(5.11) = 2n- 2c,(2n)2/3 

and scale yi according to 



(5.12) = (25,)'/' + 



Si 

For large n and an{c,y) := e-(2")'^'>'(4n)-^(2n)^/^g(2n)V3,2^2cV3^ 

(5.13) ^^^ms.^yMsum)) - ^|:f ^-^fe.>S-),(c.r.)) 
where K^°^^ is given by iS.5]) . and hence 

(5.14) Jim ( ^J^^i/g ) Pw({(5,,y,)b=i,...,.) = det[X-f*((c„y,),(c,,yO)k/,=i,...,.. 

Proof. In the case cj < ci, the formula (|5.9p still applies, but now with q 2ci(2n)^/'^. This 
latter point means that (|5.10p is now inappropriate, following |FN08[ eq. (5.40)] we use instead 
the large s expansion 



V (s - fc;)!_ 

with s — 2n, hi = Ci + 2k. A Riemann sum approximating the first integral in (|1.6p is obtained, 
establishing (j5.13p in the case Cj < ci. 

In the case Cj > c;, it follows from (|4.42p that (j5.9p remains true, but with k now ranging from 
— oo to 0, and the RHS multiplied by —1. Because the resulting Riemann integral is convergent, 
and is precisely the same as that obtained for Cj < ci except that the range is over (— cx),0), the 
second integral in (|1.6p is obtained, establishing (|5.13p for Cj > ci. □ 

5.4. Hard edge and bulk. The hard edge scaling is obtained by change of variables 
(5.16) - ^ 



2V" 

so that the mean particle density in the neighbourhood of the origin is of order unity. We seek to 
analyze the correlations with this scaling, and the species differing from 2n as specified by ()5.4p . 
Taking Yj — > oo in ths expression, with the differences between the Yj fixed, the bulk correlations 
follow as a limit of the hard edge correlations. 

Proposition 5.5. For the hard edge scaling specified by 15. J 6]) . |5.^[ ), 

(5.17) ^K{.s„y,-suyi) - 2^'-^^n^^''^^^/'K'^^^-'{{s,,y,),{si,yi)) 

where 

I 2 cos (ttwx — TT (^^)) cos (ttuj/ — TT {^^)) du, t>s 

I —2 cos (ttmx — TT {^-^)) cos (TTuy — tt (^-j^)) du, t < s 



(5.18) K^^^-''iis,x),it,y)) 
Hence 



(5.19) Jim (^^^ pi{is,,y,)},=^_r)=det[K'^^^-H{c,,Y,),{ci,Yi))],^^^^^^^ 
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Furthermore, with Sj again given by ^5.4^, the bulk scaling is obtained by the change of variables 



ttY- 

(5.20) z/i = — L + a, 0<a<2Vn 

2y/n 

and in this scaling limit 

(5.21) lim (-^) p({(s,,%)},=i,...,.)=det[A^b-d((c,,r,),(Q,ll))],-i=i,...,,. 

Proof. Consider the hard edge scahng, and suppose ci > cj. In (j5.7p we substitute for the 
Hcrmite polynomials the uniform asymptotic expansion |Sze75j 

(5.22) r(n/2 + n _^2/,^^ ^ ^^^ ^——^ _ ^ 0(n-i/2) 

i (n + 1) 

to deduce 

2 1 (sj -2fc)! 



(5.23) K{{s„y,),{si,yi))^^Y. ^- 



- 2'c(s^./2-fc)!(s,/2-fc)!- 



, n — k^^ /l^c,\\ / n — k^^ /I — c; 
X cos I n\l Yj — TT — — ^ cos ttW Y; — tt 



n 



Substituting in this (j5.4p . where still to be done, and scaling the summand by 
(5,24) 2*-«(,(n-.))-V^|^ 



(which for large n — k tends to unity) shows 

„ n 
(5.25) Ki{s„y,),isi,yi)) ^ _2^'-.„(-i+=-=.)/2 ^ 



fc=i 



(-l+c,-c,)/2 



X 



X cos -^Y, - n j j cos (^.^ --F, - n 

This is a Riemann sum to the first integral in (|5.18p , implying the leading asymptotics (|5.17p in 
this case q > Cj. 

Consider next the hard edge scaling in the case c; < cj . As in the analogous stage of the proof 
of proposition 15. 2[ it follows from (14.42p that (j5.23p remains true but with the summation now 
over (—00,0], and the RHS multiplied by —1. This gives a Riemann sum approximation to the 
second integral in (|5.18p and establishes (|5.17p in the case q < Cj. 

The result (|5.2ip can be deduced from (|5.19p by noting that the sought bulk scaling correlation 
kernel must be the limit x,y —> oo, \x — y\ fixed of K^'^'-''^ {{s , x); {t, y)). A simple rewrite of the 
resulting integrals gives (|1.7p . □ 
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